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This paper investigates the mobility of electrons scattering from the coupled system of electrons 
and longitudinal optical (LO) phonons in n-type GaAs. The Boltzmann equation is solved exactly for 
low electric fields by an iterative method, including electron-electron and electron-LO phonon scat- 
tering dynamically screened in the random-phase-approximation (RPA) . The LO phonon self-energy 
is treated in the plasmon-pole approximation. Scattering from ionized impurities screened in static 
RPA is calculated with phase-shift cross sections, and scattering from RPA screened deformation 
\ potential and piezoelectric acoustic phonons is included in the elastic approximation. The results 

ON ■ show that dynamic screening and plasmon-phonon coupling significantly modify inelastic scattering 

at low temperatures and densities. The effect on mobility is obscured by ionized impurity scattering 
in conventionally doped material, but should be important in modulation doped structures. For 
uncompensated bulk n-type GaAs, the RPA phase-shift model for electron-impurity scattering gives 
lower drift mobilities than the standard Thomas-Fermi or Born calculations which are high compared 
to experiment. Electron-electron scattering lowers the mobility further, giving improved agreement 
ON ' with experiment though discrepancies persist at high donor concentrations (n > 10 18 cm~ 3 ). When 

impurities are ignored, inelastic scattering from the coupled electron-phonon system is the strongest 
scattering mechanism at 77 K for moderate doping. This result differs from the standard model 
neglecting mode coupling and electron-electron scattering which has the acoustic modes dominant 
in this regime. 
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I. INTRODUCTION 



Good quantitative agreement between previous calculations and experiment on low-field transport coefficients of 
most high-purity polar semiconductors Q suggests that the basic electron interactions are well understood in these 
materials. This is not the case for doped polar semiconductors. A long-standing discrepancy between theoretical 
and experimental mobility values for moderately to heavily doped n-type GaAs has been noted by many authors 
P-|l3|. The fact that calculated mobilities are consistently too high, especially for moderate to heavy doping and 



low temperatures, has motivated suggestions that compensation effects 0, correlated impurity distributions [Ipl, 



modifications of screening theory for a multi-ion system jj| , or multiple scattering from ionized impurities ]l3| should 
O ■ be included to obtain accurate results. 

J-^ ' At the same time, it is known that electron-electron interactions can strongly influence electron transport in doped 
^ . semiconductors, by screening matrix elements and renormalizing phonon energies, and through electron-electron 
scattering. More accurate treatments of electron-electron interactions frequently improve agreement between theory 
and experiment. For example, screening of ionized impurity potentials in the static RPA rather than the commonly 

14| to reduce theoretical electron mobility in 
y screened electron-electron scattering gives 



used linearized Thomas-Fermi approximation (LTFA) has been shown [ 
n + silicon by more than a factor of two at 77 K. Including dynamical 

a significant reduction in drift velocity of hot electrons injected into n-GaAs devices |L5|]. Calculations including 
the coupling between plasmons and longitudinal optical (LO) phonons give shorter electron inelastic lifetimes [|L6|JlT|| 
and enhanced hot-electron energy relaxation in n-type GaAs. Previous attempts to include electron-electron 

scattering in GaAs mobility studies have been limited to static Thomas-Fermi screening of single-particle scattering 
B. So far, a treatment of inelastic scattering from the nonequilibrium coupled mode system in combination with an 
accurate model of ionized impurity scattering has not been attempted. 

In reference [ ^o| , the inelastic collision term in the electron Boltzmann equation was derived by making the Born 
approximation and the random-phase-approximation (RPA) for scattering from the coupled electron-LO phonon 
system in a doped polar semiconductor. The result is the sum of an electron-electron collision term and an electron- 
LO phonon collision term that includes plasmon-phonon mode coupling. Both interactions are dynamically screened 
by only the electronic part of the total dielectric function for the electron-phonon system. The present paper describes 
an iterative method for exactly solving the low-field Boltzmann equation including the dynamically screened inelastic 
collision terms for arbitrary electron degeneracy and spherical energy surfaces. 

For steady-state transport, electron-plasmon scattering is correctly accounted for by dynamically screening the 
electron-electron collision term pQ| . Electron-plasmon and electron-electron scattering events conserve total electron 
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momentum and energy in direct gap doped semiconductors since Umklapp processes are negligible in these materials. 
Therefore, the purely electronic scattering events cannot by themselves degrade an electrical curent, if spherical energy 
surfaces are assumed, and they merely rearrange the nonequilibrium electron momentum distribution produced by an 
applied field, thereby changing the probabilities for other scattering processes. 

In the present work, the effect of electronic scattering is accounted for by solving the Boltzmann equation for the 
nonequilibrium electron distribution with electron-electron, electron-phonon, and electron-impurity scattering treated 
simultaneously. An isotropic effective electron mass and parabolic energy bands are assumed. The phonon distribution 
is approximated by its equilibrium form, and the plasmon-pole approximation f2l] , ^2|| is used in the LO phonon self- 
energy. Numerical results are presented for electron drift mobility in uncompensated n-type GaAs, including also 
electron scattering from piezoelectric and deformation potential acoustic phonons screened in the static RPA. The 
calculations for conventionally doped material include phase-shift scattering from ionized impurity potentials screened 
in the static RPA. This method gives lower mobilities than LTFA or Born calculations of electron-impurity scattering, 
in better agreement with experiment. Results for n-type GaAs with no impurities present show that electron-electron 
scattering and mode coupling rather dramatically alter inelastic scattering at 77 K, suggesting that these effects are 
significant for modulation doped material at low temperature. 

This paper is organized as follows. Section II discusses the screening functions appropriate for the various electron 
scattering mechanisms included in this study. Section III describes the iterative method used to solve the Boltzmann 
equation. Section IV presents the numerical results with discussion emphasizing the differences obtained by using 
RPA vs LTFA screening, and by including electron-electron scattering and plasmon-LO phonon coupling. 



II. SCREENING IN A MULTICOMPONENT SYSTEM 

In reference |2p[] , it was shown that the nonequilibrium character of screened electron interactions is essential to 
obtaining the correct form of the collision terms in the Boltzmann equation, which contain nonequilibrium distribution 
functions for the system's excitations. When this is done, there remain dielectric functions which screen the interaction 
matrix elements in the collision terms. When the Boltzmann equation is linearized with respect to electric field 
strength, only terms with equilibrium dielectric functions screening the matrix elements in the collision terms are 
nonvanishing. Therefore, consider the equilibrium screening process. The interaction between two test charges in 
a solid can be expressed as the bare Coulomb interaction v q = Aire 2 /q 2 screened by the total dynamic dielectric 
function €t(q,u>) for the system, including contributions from both the conduction and lattice charges. The frequency 
dependence of the screened interaction is determined by the energy difference between the scattering electron's final 
and initial states. 

In the RPA, the electron-electron interaction is equivalent to the test charge-test charge interaction. In the same 
approximation, the total interaction v q /eT{q, u>) is separable into a purely electron-electron interaction plus an electron- 
phonon interaction including a phonon self-energy term, each of these interactions being screened by only the electronic 
part of ct- The screened electron-phonon interaction V sc —ph includes the interactions with all relevant phonons. It 
can be determined [ p3[ by subtracting from the total interaction the purely electronic interaction / e(q,uj), 

v v°° 

V sc - ph (q,w) = -1 - -2- (1) 
er e 

= Ysh. (21 

e[l-{v^ + V ph )Py { > 

where = i> 9 /eoc, e(<7,w) = 1 — v q X) P(q,Lu) is the electronic RPA dielectric function, and P denotes the polarization 
of the noninteracting electron gas in this paper. The bare electron-phonon interaction V p h is the sum over all 
phonon modes A of the product of the phonon Green's function D\ and the interaction matrix element M\ squared, 
Vph = Ylx I-^a| 2 -Da(<7,^)- This method was used in reference [Q in order to determine the screened electron-LO 
phonon interaction V sc -lo f° r the case when only LO phonons are present (Vph = Vlo)- 

When acoustic phonons are present as well, V p h = Vlo + V ac and the total screened electron-phonon interaction 
V sc -ph can be separated into V sc -lo ( as determined before) plus the screened electron interactions with acoustic 
phonons: V sc - ph = V sc - LO + V sc - ac . Solving for V sc _ oc , 

T/ / \ Vlo + Vac Vlo 

K'.c-acW.Wj - _ ( ^ + + Vac)p] £[1 _ + Vlq)p] (6) 

1 "' (4) 



(e - V LO P - V ac P)(e - V LO P) 
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Since the acoustic phonon frequencies are small for relevant q values, V sc - ac is evaluated in the static limit. Using 
Vlo{<1, 10 = 0) = (^(e^ 1 — e^o 1 ); which implies that e — VlqP = 1 — v qP/ e o- we have 



sc ac (e - v q PY[\ - e V ac P/(e - v q P)} ' v ' 

Thus, eo appears in the acoustic mode screening function. The acoustic phonon self-energy eoU ac P/(eo — v q P) is 
neglected in the present work. 

The electron-impurity interaction may also be viewed as a test charge-test charge interaction. Since the scattering 
process is assumed to be nearly elastic, the relevant screening function is er(q,^) evaluated at id = 0. In this case, 
ct should include the electronic and LO phonon parts plus a contribution from the piezoelectric phonon interaction. 
The piezoelectric contribution is assumed to be small and is neglected here so that the screening function for the 
electron-impurity interaction is eo — v q P, the same as for the acoustic modes. 



III. NUMERICAL METHOD 



In order to solve the electron Boltzmann equation for a weak electric field F, the nonequilibrium distribution 
function /(k^) may be expanded around the equilibrium solution /? = {1 — cxp[f3(Ei — to first order in F, 

/(k i ) = /° + /0(l-/P)a ; ^i, (6) 

where Xi is the cosine of the angle between F and k^, and ipi = ip(ki) is a function that is linear in F, but is otherwise 
unknown and to be determined by solving the Boltzmann equation. 

Rode jl]] devised and implemented an iterative procedure for solving the Boltzmann equation to determine /(k) to 
first order in F for the case of spherical energy surfaces, including electron-phonon and electron-impurity scattering 
only. Assuming the linear form of / (^|) and using the law of cosines in the electron Boltzmann equation (equation 
(40) of reference pp[|) gives an integral equation for the nonequilibrium electron distribution function |l]]. 

pF 8f° ■ r, 

~h&k Xk = -V^Z ar = ^ $> P a^{W(k,p)[l - /°] + lU(k,p)/°} 

p 

p 

where = /°(1 — fk)ipk, %k P is the cosine of the angle between k and p, and W(k,p) is the total differential 
scattering rate for electronic transitions from state |k) to state |p) due to scattering from ionized impurities and 
phonons. In the present work, TU(k, p) for the interaction with LO phonons is modified to include the plasmon- 
phonon coupling as described in reference |po| . Approximating the phonon distribution function by its equilibrium 
form N°(lu) — [e^ nu) — and using the plasmon-pole approximation on the phonon self-energy gives 

_ 2M 2 
yth\e(q,uJk,p)\ 2 



where M? = v q {e ^ — e 1 )huj\ /2ujTo an d the frequencies Q pi u;+, and are given by 



S$ = wi[l-e- 1 (q,0)]- 1 (10) 
i{a;l + ^±[(a;! -^) 2 + 4^K -4 )] 1/2 } . (11) 



The weight factors (w 2 . — £jp)(ui+ — lo 2 _) 1 in Im[D ] give the phonon strength in each of the hybrid uj± modes, so that 
the differential scattering rate W LO is the rate for scattering from only the phonon component of the hybrid modes. 
The screening function for the LO phonons is the electronic part of the total dielectric function, e(q, u) = l—v q K 'P(q 1 uj). 

Rode's iterative method can be generalized to include electron-electron scattering by adding the linearized version 
of the electron-electron collision term {f}^, u (equation (47) of reference 1) to the right-hand side of (0) before solving 
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for ipk- The appendix shows how this is done by again using the law of cosines and simplifying {f} e c e u- Then, the 
noncquilibrium function tp satisfies the integral equation, 



pF 8f° 



Tee(k) 



1 



1 



Tei(fc) T LO {k) 



(12) 



with the terms defined in this section. 

For elastic scattering processes, the identities k — p and W(k, p) = W(p, k) apply in equation (|^) so that the elastic 
scattering rate has the form 



w = ^E( 1 -^)^( k 'P) 



(13) 



Interactions with acoustic phonons are approximated as elastic so that rj 1 is composed of the rates for acoustic 
deformation potential, piezoelectric, and ionized impurity scattering, Tj 1 (fc) = rj p (k) + T pe 1 (fc) + (k). The 
acoustic mode scattering rates are calculated as in reference Q except the interactions are screened in the static RPA 
as described in the previous section. 



1 _ E 2 k B Tm* 

T dp (k) 4:TTh 3 Clk 3 

_ e 2 P 2 e k B Tm* 



2 k 



dq 



(l-v q P(q,0)/e o y 



T pe{k) 



eoh-'k 3 



2k 



dq 



(l-v q P(q,0)/e ) 



(14) 
(15) 



where E\ is the acoustic deformation potential, Ci is the spherically averaged elastic constant for longitudinal modes, 
and P pe is the piezoelectric coefficient. 

The electron- impurity scattering rate 1/t|^ s is calculated with the phase shifts 6i(k) determined by numerically 
solving the radial Schrodinger equation with an impurity potential screened with the static RPA total dielectric 
function eriq^O) JlJ]. The results presented below include a comparison with the Born approximation l/r i f orrl . 
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+ l>in 2 [6i(k) - 5i +1 (k)} 



1=0 



l (k) el{hkf 



N, 



2 k 



dq- 



q(l~v q P(q,0)/e ) 2 ' 



(16) 
(17) 



where Ni is the ionized impurity concentration. 

The LTFA is made by using the temperature-dependent Thomas-Fermi dielectric function 6ltfa = 1 + 
where q\ F = 47re 2 / 'e(dn/ d/i) and \i is the chemical potential. The dielectric constant e is t x for electron-electron 
or electron-LO phonon scattering, while e = eo for the screened impurity and acoustic phonon potentials. Using the 
LTFA in equation ( |l7| ) gives the familiar Brooks-Herring formula. 

The inverse lifetimes T^}(k) and T^ {k) describe the rate at which the nonequilibrium population of state |k) 
decays due to inelastic scattering with the equilibrium population of the other states. For electron-electron collisions, 



\k) 



eLhk 



dppil-fp^N ^ 



f k +P dq lm[P(q,tu k:P )} 



For scattering with the LO phonon component of the hybrid modes 

r k+k 



l\k-k~\ 



k- P \ q 3 Hq,^k,pW 

dql(q,uj x ) N°(uj x ) + 1 - f k -]®(Ek - Tw X ) 



(18) 



k+kl 



dql(q,u x ) N°(uj x ) + f° 



fc-fc+| 



(19) 



where 



x/k 2 ±2m*Lu x /h, 



e uj LO m 
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Since the mobility is measured on a system with all of the electrons out of equilibrium, the linear integral equation 
( |l2l ) also contains the functional tWi^A] describing the rate of change of the equilibrium occupation of the state |k) 
due to inelastic scattering with the nonequilibrium population of the other states. This term is composed of electron- 
electron and electron-LO phonon parts: Vi ne i = v CP + vlo- The contribution from electron-electron scattering is 



-4e 4 / fc ° ' v 



r k+p fin 

dp (1-/°) / 4\e(q,u; k , p ) \~ 

J \k—p\ H 



e^hk 2 

x \ M k2 +P 2 - 1 2 W°(uk, P ) + l}Im[P(q,LJ kiP )} 



; (k -p +q ) 



2tt% q 



2 n2 



\*- 



ds </> s / s °(l - /° + )-z+ ds ^ s /°- (1 - / s °) 

J\z+\ 



(20) 



where z ± = (k 2 — p 2 ± q 2 )/2q and s ± = \/ s 2 ± (k 2 — p 2 ). The contribution from the LO phonon strength in the 
hybrid modes is 



k+k. 



dql(q,uj\)g(k x ) 



(k 2 + (k- x ) 2 -q 2 ) 



\k-k~\ 
k+k+ 

k—ki 



N°(Lu x ) + f° e(E k -Hu; x ) 



dql(q,ux)g(kx) 



(k 2 + (ktf - q 2 ) 



N°(co x ) + l-f k 



(21) 



Plasmon-phonon coupling may be neglected by setting w+ equal to u>lo an d UJ- equal to Q p in equations ( |l9|) and 



IV. RESULTS AND DISCUSSION 

The numerical results for the drift mobility were obtained by solving the iterative equation ( |l2| ) with (^) for the 
linearized nonequilibrium electron distribution in order to determine the average velocity per unit electric field. The 
effective electron mass value m* — 0.07m e was used throughout the calculations. The other material parameters 
used were taken from reference p5[ . Results for uncompensated n-type GaAs are presented for conventionally doped 
material (figures 1 through 4) and for a doped system with no impurities present, an idealized model of modulation 
doped material (figures 5 through 7). 



A. RPA vs LTFA screening 



The LTFA is an approximation to the static RPA dielectric function which does not account for the polarization 
of the screening electrons by the colliding electron, but is equivalent in the limit of small momentum transfer, q. 
Figures 1 and 2 show the carrier concentration dependence of electron mobility in uncompensated n-type GaAs at 
300 and 77 K, calculated with statically screened electron-phonon and electron-impurity interactions. Scattering from 
LO phonons as well as the acoustic piezoelectric and deformation potential modes was included, but electron-electron 
scattering and plasmon-phonon coupling were neglected. Each figure shows calculations made by using either the 
phase-shift (|l6| ) or Born (|l7]) electron-impurity scattering rate, and either the temperature-dependent static RPA or 
LTFA dielectric function in equations ( fbif - 17] ) . The standard LTFA Born (Brooks-Herring) calculation with unscreened 
phonons is also shown. Compared to the case of n-type Si fll^| , a relatively small correction to mobility in n-type 
GaAs is obtained by using RPA screening rather than LTFA. Except for the unscreened phonon case at 300 K, the 
combination of RPA screening of impurities and phonons with the phase-shift electon-impurity scattering rate gives 
the lowest mobilities, while the LTFA Born calculations give values that are higher by at most 16% at 300 K and 
41% at 77 K. A cancellation of errors made by using the LTFA Born electron-impurity interaction and unscreened 
electron-phonon interactions at 300 K makes this result similar in magnitude to the RPA Born curve. Dynamic RPA 
screening of LO phonons changes the mobility very little from the static RPA case when electron-impurity scattering 
is included. 
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B. Electron-electron scattering 



The effect of electron-electron scattering on mobility in doped semiconductors is limited by three factors. 1) The 
screened interaction between two electrons is considerably weaker than the Coulomb interaction. 2) When Umklapp 
electronic scattering processes are negligible and energy surfaces are close to spherical, electron-electon scattering 
conserves total electron momentum and produces only an indirect effect on other scattering processes by redistributing 
the electron momenta. 3) The Pauli principle restricts the fraction of electrons that can participate in energy- 
conserving electron-electron scattering processes to a number that vanishes in the degenerate limit. Nevertheless, 
electron-electron scattering was shown \2& to have an important effect on transport in doped semiconductors away 
from the degenerate limit. Chattopadhyay Q studied its effect on mobility in n-type GaAs by using a variational 
solution of the Boltzmann equation including electron-electron scattering statically screened in LTFA with eo- He 
found that, for n — 10 16 cm -3 , electron-electron scattering reduces the mobility by about 10% at 80 K. The reduction 
increases to 20% when rather than eo is used in the LTFA screened electron-electron interaction, as discussed 
in reference pC( |. Figures 3 and 4 show the drift mobility calculated by iteratively solving the Boltzmann equation 
at 300 and 77 K, including electron-electron scattering in addition to the electron-impurity and electron-phonon 
mechanisms already discussed. The Pauli principle restriction is clearly seen from the vanishing of the electron- 
electron contribution at higher concentrations. Also shown in figures 3 and 4 are experimental data |6|,p|p7|] for the 
Hall mobility of n-type GaAs. No attempt has been made in the calculations to include the Hall factor relating the 
drift to Hall mobilities. The Hall factor is generally larger than unity so that the Hall mobility is larger than the drift 
mobility, but its magnitude depends on doping density and temperature, approaching unity in the degenerate limit. 
Despite the reduced mobility values compared to previous calculations (see e.g. [|l|,||J^J^] ) , a disagreement between 
theory and experiment persists for carrier concentrations above n = 10 18 cm -3 at both temperatures considered. 



C. Coupled plasmon-phonon modes 

The uncoupled plasmon temperature Tiuip/kB varies as n 1 / 2 ; for n-type GaAs, it ranges from 50 K to 1594 K for 
electron concentrations of 10 16 cm~ 3 to 10 19 cm~ 3 . The plasmon temperature crosses the LO phonon temperature 
of 419 K at ~ 7 x 10 17 cm -3 . In this concentration region, the electronic and lattice excitations hybridize to form 
the normal modes of mixed electron-phonon character. The other energy relevant to electron mobility is the Fermi 
energy Ef, which varies as n 2 / 3 . For n-type GaAs, EfI^b ranges from 29 K to 2930 K for n — 10 16 cm -3 to 10 19 
cm~ 3 and crosses the LO phonon temperature at ~ 6 x 10 17 cm~ 3 . 

In the plasmon-pole model, the phonon spectral function — 7r _1 Im[_D(g, oj)] has two delta-function peaks (see equa- 
tion (|9|)) with frequencies tu + and given by (pd[). As discussed in reference mobilities at low densities and 
temperatures are expected to be reduced when plasmon-LO phonon coupling is included, because of increased low- 
energy electron-phonon scattering due to the u>- mode. Figure 5 shows that this is in fact the case. The mobilities in 
this figure were calculated by neglecting all scattering mechanisms except electron-LO phonon scattering. For a wide 
range of densities at 77 K, the LO phonon limited mobility including mode coupling and dynamic RPA screening 
(solid curve) is significantly lower than the mobilites calculated by neglecting these effects. The dip observed below 
n = 10 18 cm~ 3 in all the figure 5 mobility curves is related to Ep crossing the LO phonon emission threshold. An 
interesting though minor effect appears in figure 3, where the mobility including mode coupling is actually higher at 
low doping levels and 300 K than the mobility neglecting the coupling, when electron-impurity is included. 

Since mobility is limited primarily by scattering from ionized impurities at low temperatures in conventionally doped 
GaAs, the mode-coupling effect is obscured in this case. Nevertheless, it could be important for mobility in modulation 
doped structures where the effects of impurity scattering are greatly reduced because of the large separation between 
dopants and carriers. Figure 6 plots the electron drift mobility as a function of concentration at 77 K when electron- 
impurity scattering is neglected. The figure shows that scattering from the acoustic (piezoelectric and deformation 
potential) modes is stronger than scattering from LO phonons when electron-electron scattering is neglected, even if 
mode coupling is included. However, in the moderate doping regime, inelastic scattering is stronger than scattering 
from acoustic phonons when both electron-electron scattering and mode coupling are included. This result is notable 
especially because the value used in the calculation for the acoustic deformational potential E\ — 12 eV is higher 
than some others which appear in the literature. Figure 7 plots the temperature dependence of the mobility from 50 
to 150 K for n = 10 17 cm -3 when electron-impurity scattering is neglected. The figure shows that electron-electron 
scattering and mode coupling effects should be significant even at the lower temperatures in situations where impurity 
effects are minimal. 
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V. CONCLUSIONS 



The mobility calculations presented in this paper show that a more accurate treatment of electron interactions 
in n-type GaAs yields results that differ from simple models, especially at moderate doping levels in conventionally 
doped material and for low densities and temperatures in modulation-doped material. For conventional n-type GaAs, 
RPA screening of phonons and impurities combined with phase-shift electron-impurity scattering rates give lower 
drift mobilities than LTFA screening or Born calculations. A convergent iterative method has been developed for 
solving the Boltzmann equation including electron-electron scattering. Including this scattering mechanism further 
reduces the mobility except at high doping levels, so that the full calculation gives better agreement with experimental 
Hall mobilities than previous results. In order to say definitively that the long-standing discrepancy between theory 
and experiment is resolved in the middle concentration regime, the Hall mobility should be calculated with the 
present model. At the highest concentrations, the theoretical mobility for uncompensated n-type GaAs is larger than 
experimental values. 

The effects of dynamic screening, mode coupling, and electron-electron scattering are not very important for mo- 
bility in conventionally doped GaAs where scattering from ionized impurities dominates. However, the calculations 
neglecting impurities suggest that these effects could significantly affect mobilities in modulation-doped structures. 
Comparison with experiment would require a more realistic treatment including the effects of confinement on phonon 
spectra, the two-dimensional nature of the electron gas in these structures, as well as the effects of impurities, alloy 
disorder, and interface roughness. Also, the effect of Landau damping of the hybrid modes, which is neglected in 
the plasmon-pole model for the phonon self-energy, should be investigated. Good candidates for study would be 
ALcGai^As, since the Frohlich coupling increases with x pSj , or the more strongly polar materials CdTe, ZnSe, or 
ZnTe. Finally, the fact that inelastic scattering is stronger compared to acoustic phonon scattering than the standard 
model predicts for doped material at low temperature suggests that the inelastic component could influence deter- 
minations of the deformation potential constant from mobility studies, and motivates a reevaluation of Ei in n-type 
GaAs at low temperatures. 
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APPENDIX: ITERATIVE METHOD FOR ELECTRON-ELECTRON SCATTERING 



Rode's method Q is generalized by inserting the linearized form of / (^|) into the electron-electron collision term 
(equation (47) of reference [51) and using the identity 

5{E X + E 2 -E 3 - E 4 ) - - /°) - (1 - /°)(1 - / 2 )/ 3 °/ 4 ] = 0, (Al) 

giving the linearized collision term, 

{/ (ki)}£« = 7^3 E l M 3,i l^(Si + ^2 - E 3 - £? 4 )5(ki + k 2 - k 3 - k 4 ) 

2,3,4 

X /f/2 (1 - /aX 1 - /4 )NlV>l + ^2 - - X 4 ^ 4 } • 

This symmetric form allows a simple application of the law of cosines to eliminate x$ for i = 2 to 4 in favor of Xn, 
the cosine of the angle between ki and k;. 



{/ (ki)}S« = E l M 3,i I 2 ^(^i + E 2 - E 3 - £7 4 )*(k! + k 2 - k 3 - k 4 ) 



X A°/2°(l - / 3 °)(1 - /4°)[^1 + ^12^2 " X 13 ^3 - X U M 



(A2) 



Now use the delta function for momentum conservation to sum over k 4 and introduce the variable q = k 4 — k 3 = k 4 — k 2 . 
The k 2 integration of the ip\ and ip 3 terms can be done analytically. Making the variable changes ki — ► k, k 2 — > 
s, k 3 — > p — k — q, and k 4 — > s + q, we have 



(kl !■':// = ^ 1 Mkp |2 ^~ x ^i N °M + l]MP(<W P )]/ fc (l - f p ) 



~ 2 I M fe , p | 2 x fes 5 (P fc - P p + ?^ )S+(? )/°/°(l - - /?+,) 



where s was changed to s — q in the last term, and 



Im[P(g )W )] = — w 
2nn pq 



In 



1 + exp[-(3(E+ - im)] 



1 + exp[-0(E~ - M) 



(A3) 



(A4) 



with E ± = (E q ± hu) 2 /AE q . 

The symmetry of the problem is best used by choosing bipolar coordinates j2j|] with the z axis along q and the xz 
plane containing k. The p integration variables can be chosen to include the momentum transfer q, 



k+p 



d p = 2n / dpp / dq — , 

Jo J\k-p\ k 

the factor of 2ir coming from integration of the azimuthal angle around k. The s integration has the form J d 3 s — 
J dcj) J dx sq J ds s 2 where <p is the angle between the (k, — p) and (— s, k 4 ) plane. The <f> integral is simply done by 
again applying the law of cosines, 



2 IT 



' ^ks ^7T X sq Xf~ q . 



Then, using the energy delta functions to do the x sq integrals, the linearized electron-electron collision integral reduces 
to 

"k+p 



if (k)} 



-x k f? r°° 



dpp(l - p I dq I M Kp Y 

'k—p\ 



U M2^) 3 Jo 
x\ q (ijj k - x kp ?pp)[N (uj k ,p) + l]Im[P(g 



+ 



2ttTi 



z~ / ds - / s °+ ) - z+ / ds ^ s / s °_ (1 - / s °) 



(A5) 



8 



where z ± — (k 2 — p 2 ± q 2 )/2q and s^ 1 — yjs 2 ± (k 2 — p 2 ). 

Adding (AE) to the right-hand side of (Q) and solving for tpk gives equation 
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FIG. 1. Comparison of RPA vs LTFA screening and phase-shift vs Born electron- impurity cross sections at 300 K for n-type 
GaAs. Static screening of LO, piezoelectric, and acoustic deformation potential phonons as well as ionized impurities is included. 
The RPA phase-shift drift mobility (solid curve) is shown with mobilities in the following approximations: RPA Born (dotted 
curve), LTFA phase-shift (long-dashed curve), LTFA Born (short-dashed curve), and LTFA Born with unscreened phonons 
(chain-dashed curve). Electron-electron scattering and plasmon-phonon mode coupling has been neglected. 

FIG. 2. Comparison of RPA vs LTFA screening and phase-shift vs Born electron-impurity cross sections at 77 K for n-type 
GaAs. Calculated mobilities are represented as in figure 1. 

FIG. 3. Effect of electron-electron scattering on electron drift mobility in n-type GaAs at 300 K. The dotted and dashed 
curves are the RPA phase-shift calculation neglecting and including, respectively, statically screened electron-electron scatter- 
ing. The solid curve includes dynamically screened electron-electron and electron-LO phonon scattering with mode coupling. 
Experimental Hall mobilities are indicated by O ( re f- [9]) an< i '-' ( re f- [27]). 

FIG. 4. Effect of electron-electron scattering on electron drift mobility in n-type GaAs at 77 K. Calculated mobilities are 
represented as in figure 4. Experimental Hall mobilities are indicated by O ( re f- [9]) an d A (ref. [12]). 







FIG. 5. Effects of screening and mode coupling on LO phonon limited electron mobility in n-type GaAs at 77 K. The 
dynamically screened coupled mode mobility calculation (solid curve) is compared to mobilities determined by scattering from 
uncoupled LO phonons screened with LTFA( long- dashed curve), static RPA (short-dashed curve), dynamic RPA (dotted curve), 
and unscreened (chain-dashed curve) . All other scattering mechanisms are neglected. 

FIG. 6. Comparison of mobilities limited by inelastic and acoustic phonon scattering in n-type GaAs at 77 K. The mobility 
including dynamically screened electron-electron and electron-LO phonon scattering including mode coupling (chain-dashed 
curve) is compared to the same calculation except neglecting electron-electron scattering (dashed curve), and to the mobility 
including only scattering from acoustic (piezoelectric and deformation potential) phonons screened with static RPA (dotted 
curve). The solid curve shows the calculated mobility including all of the mechanisms. Electron- impurity scattering is not 
included. 

FIG. 7. Temperature dependence of mobility in n-type GaAs for n = 10 17 cm -3 . The mobility including dynamically screened 
electron-electron and electron-LO phonon scattering including mode coupling (solid curve) is compared to the results with static 
RPA screening of electron-electron scattering (dotted curve), electron-electron scattering neglected (dashed curve), and static 
RPA screening of LO phonons without mode coupling (chain-dashed curve) . Electron scattering from acoustic phonons screened 
in static RPA is included and electron-impurity scattering is neglected in all cases. 
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